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Abstract: Small-x evolution in QCD is conveniently described by Mueller's dipole model 
which, however, does not include saturation effects in a way consistent with boost invari- 
ance. In this paper we first show that the recently studied zero and one dimensional toy 
models exhibiting saturation and explicit boost invariance can be interpreted in terms pos- 
itive definite k ^ k + \ dipole vertices. Such k ^ k + \ vertices can in the full model be 
generated by combining the usual dipole splitting with k — \ simultaneous dipole swings. 
We show that, for a system consisting of N dipoles, one needs to combine the dipole split- 
ting with at most — 1 simultaneous swings in order to generate all colour correlations 
induced by the multiple dipole interactions. 
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1. Introduction 

The small-x region in QCD can be described by the well known, linear, BFKL equation 
which predicts a power like growth in x for the gluon density. Such a fast growth is 
problematic since it breaks the unitarity bound at high energies. As the gluon density 
becomes large, non-linear effects cannot be ignored and it was early suggested that effects 
from parton saturation should tame the growth of the gluon density, in accordance with 
unitarity [1]. 

Starting the evolution from a colour singlet quark-antiquark pair, a colour dipole, 
Mueller [2,3] formulated a dipole model in transverse coordinate space which reproduces 
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the BFKL equation to leading order. The transverse coordinate formulation also allows 
one to go beyond the BFKL equation since it is here easier to take into account multiple 
interactions. This is so because the transverse coordinates of the partons are frozen during 
the evolution, and it is therefore rather easy to sum the multiple scattering series in an 
eikonal approximation. This was exploited by Mueller, who was thus able to obtain a 
unitarised formula for the scattering amplitude. 

Within the dipole formalism, Balitsky [4] derived an infinite hierarchy of equations for 
the dipole scattering amplitudes. Kovchegov [5] derived a closed equation for the amplitude 
using a mean field approximation, and this equation is referred to as the Balitsky-Kovchegov 
(BK) equation. The same hierarchy of equations also follows from the JIMWLK equation 
[6-9] which is the master equation of the Color Glass Condensate (CGC) formalism [10] and 
describes the non-linear evolution of dense hadronic matter in the presence of saturation 
effects. 

In Mueller's model the multiple dipole interactions correspond to multiple pomeron 
exchange, and in the Lorentz frame where the collision is studied these multiple interactions 
lead to the formation of pomeron loops. However, these loops cannot be formed during 
the evolution of the dipole cascade since this evolution is linear. Thus only those loops 
which are cut in the specific Lorentz frame used for the calculation are accounted for, 
while none of the remaining loops is included. This implies that the model is not frame 
independent. To minimize the error, the optimal frame to use is the one where the colliding 
dipole cascades are of the same density, since multiple scatterings then become important 
at rapidities where one may still neglect saturation effects in the evolving dipole cascades. 

In order to obtain a frame independent formalism it is necessary to include saturation 
effects also in the evolution of each dipole cascade. There have been various attempts to 
include such saturation effects in a consistent manner, but no explicitely frame independent 
formalism has yet been presented. 

To gain insight and possible hints towards a solution, a simple 1+0 dimensional (rapid- 
ity constitutes the only dimension) toy model in which transverse coordinates are neglected 
was constructed in [3,11]. This model has been further studied in [12,13], and the resulting 
frame-independent evolution can be interpreted as a coherent emission of new "dipoles" 
from the multi-" dipole" state. 

In [14] we developed a dipole cascade model for DIS and pp collisions based on Mueller's 
model but also including effects of energy-momentum conservation. This model was fur- 
ther extended in [15] to include saturation effects in the cascade evolution through the so 
called dipole swing mechanism [15,16], which gives an additional 2 — > 2 transition during 
the evolution. Monte Carlo (MC) simulations show that the evolution is almost frame 
independent, and the model results are in good agreement with inclusive and diffractive 
data from HERA and the Tevatron [17]. 

In this paper we will first show that the explicitely frame independent evolution in the 
toy model mentioned above, and also in its 1-1-1 dimensional generalization in [18], can be 
given a probabilistic interpretation in terms of positive definite k ^ k + 1 dipole vertices. 
Such k ^ k + \ transitions can in the full model be generated by combining the dipole 
splitting with the dipole swing. 
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In case each individual dipole is restricted to single scattering only, we show that one 
needs to combine the dipole splitting with at most one swing at a time in order to generate 
the necessary colour correlations. As remarked in [12], the toy model evolutions mentioned 
above describe the multiple scatterings of individual dipoles, and we will here show that 
one can in the real model generate the correlations induced by the dipole scatterings by 
combining a splitting with several simultaneous swings. For a system consisting of N 
spatially uncorrelated dipoles, it is easy to see that one needs to combine the dipole splitting 
with at most — 1 simultaneous swings. In a process where one splitting is combined with 
k — 1 swings, k dipoles are replaced hy k + 1 dipoles, thus giving a A: — > A: + 1 transition. 

However, starting the evolution from a single qq pair, one obtains dipoles which are 
connected in chains, and in this case not all swings are allowed. It is here important to 
keep track of the correct topology of the dipole state. While this is never a problem in the 
original formulation which only contains the 1 — > 2 splitting, it is here very important to 
avoid the formation of unphysical states. Although it has been checked for a large number 
of formal proof that it is always enough with A^ — 1 swings is not available, and 

the result is a conjecture. 

The paper is organized as follows. In the next section we shortly review Mueller's 
dipole model and the question of frame independence in the evolution. In section ^ we will 
review the toy models in zero and one transverse dimensions formulated in [11, 18], and 
show how these can be interpreted in terms of the k ^ k + 1 transitions mentioned above. 
Then in section ^ we consider the evolution in the full model and argue that also in this 
case the correct evolution can be formulated in terms of — > A; + 1 vertices. In section ^ 
we go on to study the colour topology of the evolution, and we show how one can generate 
the needed colour correlations using the dipole swing. Finally, in section ^, we present our 
conclusions. 

2. Approaches Towards a Frame Independent Formalism 

In Mueller's model [2,3] for onium-onium scattering a colour dipole formed by a colour 
charge at transverse coordinate x and an anti-charge at y can split into two dipoles by 
emitting a soft gluon at position z with the following probability 



= " jj2„ jx - y) 
dY 27r 



'^'^7:-^^327^^^ = '^'^^(^'y'^)' (2-1) 



a = ^ and Y = ln—. (2.2) 
vr yiBj 

We refer to A4 as the dipole kernel and to Y as the rapidity, which here acts as the time 
variable in which the evolution proceeds. A dipole cascade is then formed when each 
dipole splits repeatedly. When two such cascades collide, a right-moving dipole (xj, yi) and 
a left-moving dipole {xj,yj) interact with a probability 



fij — f{xiiyi\xj,yj) — Q 



log 



(xi - Xjf{yi - yjf 



(2.3) 
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All dipole interactions axe assumed to be independent, and the S-matrix element is given 
by 5 = exp(-^.^./ij). 

In this formalism, saturation effects occur only due to multiple scatterings while the 
evolution of the dipole cascade satisfies the usual BFKL equation. This implies that the 
formalism is not frame independent, and in order to obtain a frame independent formalism 
such saturation effects must be properly included in the cascade evolution. 

Different approaches have been proposed to obtain this. It was noted that the CGC 
formalism is not complete in the sense that it does not contain any gluon splittings, or 
gluon number fluctuations. This problem comes from the fact that the Balitsky-JIMWLK 
(B- JIMWLK) equations^ , which can schematically be written as 

dY{T^) =M® {{T^) - (r^+^)}, (2.4) 

only couples the /c— dipole scattering amplitude to the {k + n)-dipole amplitudes with 
n = 0, 1, ... . From the view of target^ evolution, this means that one includes all gluon 
merging diagrams, while the gluon splitting diagrams are absent. Gluon splittings are 
equivalent to dipole splittings, and the dipole model has been used to add fluctuations into 
the formalism. The modifled B-JIMWLK equations then read 

av-(r*) =M® {(T*) - (T'^+i)} + {T^-^), (2.5) 

where .J^ is a kernel representing the fluctuations in the target. 

Viewed in the opposite direction, the fluctuation effects in the target correspond to 
saturation effects in the projectile evolution. Given the form of the modified B-JIMWLK 
equations above, it may seem natural to include a 2 — * 1 vertex in the dipole evolution. 
Such an interpretation has, however, the drawback that it cannot be interpreted as a 
classical evolution process since the 2 — > 1 vertex is not positive definite, as was shown by 
lancu et al. [19]. In [16] Kozlov et al. calculated directly the 4 — > 2 gluon merging vertex 
within the dipole language, and were led to the conclusion that one should include a 2 ^ 4 
dipole vertex which is composed of a splitting and a swing. 

In [15] we included the 2 — > 2 dipole "swing", in addition to the 1^2 dipole splitting, 
in the evolution. If we have two dipoles {xi,yi) and {xj,yj), the swing will replace them 
by {xi,yj) and (xj,yi). The dipole swing can be interpreted in two ways. First, as a 
way to approximate colour quadrupoles as two independent dipoles formed by the closest 
charge-anti-charge pairs, in which case the swing is naturally suppressed by N"^. Secondly, 
we may view it as the result of a gluon exchange between the dipoles, which results in a 
change in the colour fiow. In this case the swing would be proportional to a^, which again, 
compared to a, is formally suppressed by A^^- 

The dipole swing is related to the pomeron interactions studied by Bartels and Ryskin 
[20,21]. Here a pomeron is interpreted as two gluons in a colour singlet state. In a four 
gluon system with two singlet pairs, gluon exchange can give a transformation ( a "switch" 

'^Throughout this paper, we will only consider the large- iVc version of these equations. 
^In the CGC approach, it is usually assumed that the target is a dense hadron while the projectile is an 
elementary colour dipole. 
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or "swing") (12) + (34) (13) + (24), where a parenthesis denotes two gluons in a singlet 
state. 

We note that the swing is here not a vertex in the same sense as the sphtting process 
since, unlike the splitting, the swing is not proportional to dY but rather happens instan- 
taneously. In the MC implementation we assign a "colour" to each dipole and two dipoles 
are allowed to swing if their colour indices match^. The swing is then determined by the 
weight 



{xi - y2y{x2 - yif 



P(swmg) oc ^. _^ (2-6) 



Here the two initial dipoles are determined by the coordinates (xi,yi) and {x2,y2), and 
the new by {xi,y2) and {x2,yi)- This form favours the formation of small dipoles. It also 
preserves one of the results in Mueller's original formulation, namely that the total weight 
for a dipole chain is given by the product Yii -^i where rj is the size of dipole i and the 
product runs over all "remaining" dipoles in the cascade. 

In our formalism the total number of dipoles does not decrease. For each event, many 
of the dipoles will not interact and these have to be considered as being virtual. In this 
case saturation effects do not have to decrease the total number of dipoles but rather only 
the number of interacting, or "real", dipoles. The dipole swing has this property since it 
is more likely that two dipoles are replaced by two smaller dipoles, as can easily be seen 
from (|2.6| ), and smaller dipoles have smaller interaction probabilities. Thus the number 
of interacting dipoles will actually decrease, and in pomeron language this means that the 
swing generates pomeron mergings. 

An essential feature of our formalism is that the dipoles in the cascade form connected 
chains^, rather than a collection of uncorrelated dipoles, as in a reaction-diffusion type of 
formalism. A dipole chain cannot end in a gluon, and it is not possible to remove a dipole 
without reconnecting its neighbors. A generic 2^1 vertex is therefore not possible in 
this formalism^, and the dipole swing gives the simplest process from which one can form 
closed chains during the evolution. 



3. The Toy Models 

3.1 The l-|-0 dimensional toy model 

In this section we will review the toy model which was studied in detail in [12] (see also [13]). 
This model was first presented by Mueller [3] and it is interesting since, besides having some 
structural aspects in common with the dipole model, it offers analytical solutions which 
have been very hard to obtain for the full model. 

^The number of effective colours is which is the number of possible colour configurations for a given 
colour-anti-colour pair. 

*In Mueller's original formulation this is not relevant since there the dipoles evolve truly independently. 
However, in our implementation of energy-momentum conservation, neighboring dipoles affect each other 
and it is then relevant that the cascade is formulated as a dipole chain. 

^The only allowed merging process is in case two neighboring dipoles merge. 
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The model is defined such that at any rapidity Y the system is specified only by the 
number of dipoles. The probability to find the system in the n-dipole state at time Y is 
denoted by Pn{Y). Here transverse coordinates are completely neglected and Y defines the 
only coordinate in the model. 

If we let 7i denote the Hamiltonian of the system we have 

'Hnm = {m[H\n) = TZ{n){Sm,n+l - Sm,n)- (3.1) 

TZ{n) is so far an unspecified function, which determines the splitting rate of the n dipole 
state. The probability Pn(Y) evolves according to 

dy Pn — 'HmnPm- ('^•^) 

With Snm we denote the S'-matrix for the scattering of two dipole states of n and m 
dipoles respectively. If we assume that each dipole scatters independently with a probability 
r, we have Snm = (1 ~ t)"^'"^. We here assume that r << 1 and if this is not the case, 
one should replace 1 — r by exp(— r). The physical 5— matrix, S, is obtained by taking an 
average over all possible events. Using matrix notation we have 

S{Yi + Y2) = p^{Yi)Sq{Y2). (3.3) 

Here p^{Yi) = {Pi{Yi), ^2(^1)1 • • • ) is the row vector of the configuration probabilities for 
the right moving onium (evolved up to Yi) while q denotes the column vector of configu- 
ration probabilities for the left moving onium (evolved up to 12)- 

In ( |3.3D we have anticipated that S depends only on the total rapidity interval Yi + 127 
which defines boost invariance. This implies that we have {dyi — dY2)S = 0, and requiring 
this in ( p.3| ) one obtains 



p^{Y^)nSq{Y2) -p^{Y,)Sn^q{Y2) = 0. (3.4) 
A sufficient condition for a solution is to require that 

ns = sn^ (3.5) 

which means that TiS is symmetric (since S is symmetric). It is now easily seen that 
condition ( |3.5| ) requires TZ{n) in ( |3.lD to be given by 

7^(n) = c(1-(1-t)") (3.6) 

where c is a constant, c = TZ{1)/t. By rescaling Y , we might as well assume c to equal 1. 

3.2 Stochastic evolution with A; — > /c + 1 vertices 
We note that (|3.6|) can be rewritten as 



(3.7) 
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This suggests that we can interpret the evolution in terms of A; — > /c + 1 transitions with 
weights (— 1)'^t'^'. However, the alternating signs implies that one cannot interpret these 
vertices in a probabilistic formulation. 

We will now show that one can nevertheless interpret the evolution in terms of positive 
definite A; — > /c + 1 vertices. Thus the evolution can still be formulated as a stochastic 
process, but we will also see that the probabilistic interpretation of the evolution implies 
that it cannot be reduced into a formalism which describes a system of incoherent particles. 

Assume we have a system of n particles X, which we also refer to as dipoles, satisfying 
the following rules. Each isolated X can emit a new X with a probability per unit time 
(rapidity) given by r, i.e. we have a reaction X ^ X + X which occurs with probability 
r > 0. In addition to this, k isolated dipoles can undergo a transition kX — > {k + l)X with 
probability . 

The evolution of the probabilities PniX) for these X then satisfies 

n n— 1 

dyPniY) = -Y.T^tlk^,Pn{Y) + Y,nt-ll,Pn-i{Y) (3.8) 

fc=l k=l 

(n) 

where TZ\_^f^j^i is the splitting rate for the process where k dipoles are replaced by A; + 1 
dipoles in a state containing n dipoles. 

Now, the splitting rate T^^"^^,,,^ is not simply given by (^t^ as one could expect 
naively, but it is instead given by 

<l+l = 0-'^(l--^-^ (3-9) 

since for each A; ^ A; + 1 we must also multiply with the probability that no more than k 
dipoles were involved in the emission of the new dipole. Obviously, k must run from 1 to 
n, and summing all contributions we obtain the total splitting rate as 

= E ( J^'(l - ^)""' = (1 - (1 - ^T) = (3.10) 

k=l k=l ^ ^ 

where TZ{n) was defined in ( |3.6D . The positive definite transition rates in (|3.9| ) thus give a 
boost invariant evolution as before. Note also that in this case the A; ^ A; + 1 transitions 
are not universal since they depend on n, unlike the rates in ( |3.7[ ). The A; ^ A; + 1 splitting 
vertex therefore not only depends on the state of the k emitters, but it does also depend 
on the rest of the dipoles in the cascade. We will in the forthcoming sections see that the 
dipole swings give a very similar evolution in the full model. 

We also note that the A: ^ A; + 1 transitions can be made manifest by writing down 
the Hamiltonian 

oo ^ oo oo 

^ = E ^(^"'^'"^ - n (1 - ^) E (3-11) 

k=l ' l=k+l m=l 

where and a are dipole creation and annihilation operators respectively, and iV" = a^a. 
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3.3 Evolution equations 

In this section we will first show that the evolution equations for the scattering amplitudes 
derived in [12] are described exactly by the /c — > A: + 1 transitions in ( |3.9| ). We will then go 
on to point out that there is a fundamental structural reason for the fact that the attempts 
to interpret the full model evolution given in ( p.5D in a probabilistic manner have run into 
problems. We will see that it is not possible to interpret this equation using a probabilistic 
2^1 vertex even in the toy model. However, we note that it is not necessary to include 
a 2 ^ 1 vertex to obtain saturation. In fact any 2 — > n vertex where only one of the n 
dipoles interact also corresponds to a 2 — > 1 transition, and this will be discussed more 
later. 

First we write the S-matrix given in ( |3.3D as 

SiYi+Y2)=p^iYi)sgiY2), Sg = Sq (3.12) 

where the nth component of the vector Sq{Y2) is the S-matrix element of a projectile, 
evolved up to Yi, made up from n dipoles scattering against a generic target, which is 
evolved up to Y2. It is then easy to see that {Sq)n = (s") satisfies the following evolution 
equation [12] 

ay(s")=7^(n){(s"+l)-(s-)} (3.13) 

Using the relation s = 1 — t, where t denotes the scattering amplitude, one can similarly 
derive the equations obeyed by (i"). Since r is assumed to be small, one can expand TZ{n) 
in each equation and drop contributions which are negligible in all regimes (dilute and 
dense systems). Doing this, the authors in [12] arrived at the following evolution equations 
for {t), (t^) and (t^), 

dy{t) = {t) - {t'), (3.14) 
dY{t') = 2{{t') - {t^)) + T{t{l - tf), (3.15) 
= 3((t3) - (t^)) + 3T(t2(l - tf) + T^{t{l - tf). (3.16) 

If one neglects all terms proportional to r, then it is seen that the resulting hierarchy 
corresponds to the large Nc version of the B-JIMWLK hierarchy. 



Let us now see how equations ( 3.14| )-( ^.16 ) arise from the transition rates TZ^"'^ 



The evolution of the S'-matrix elements are given by 

dy{s) = n^^l2{{^') - (.)}, (3.17) 

dy{s') = (TZfi, + n^^l,){{^') - is')}, (3.18) 

dris') = {nfl, + nfl, + nflMs') - {s')}. (3.19) 

It is then straightforward to obtain the following equations for the scattering amplitudes 

dY{t) = 7^£2{(^) - {f)} = it) - {t'), (3.20) 
dvit') = 2n'lU{t) - it')) - n^^Utii - tf) - n^^Uta - tf) 

= 2{{t') - it')) + 2r(t(l - t)') - r{t{l - t)'). (3.21) 
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We indeed see that the first equation is equal to (|3.14| ) and that the second equation is 
equal to ( |3.15| ). It is also straightforward to see that the equation for (t^) agrees with 

4m- 



Next, we comment on the structure of the equation given in (p.5|). Assume we wish to 
view the process as a stochastic evolution with a 1 — > 2 splitting vertex, /i^2; and a 2 — > 1 
merging vertex, A;2-+i- We here assume the total splitting rate to be the incoherent sum of 
the individual splitting rates. For the evolution of the 2— dipole state we have 

dvis') = A%{{s') - is')} + - is')} (3.22) 

which gives 

dy{t') = 2fi%{{t) - it')} + k^il,{t{l - t)) - fi%{t{l - tf) 

= (2/a. - f[% + k^^.m + i-2fi% + 2fi% - kfl,){e) - /S,(t3) 

= kfl,{t) + {fi% - kfl,){t^) - fl%{t'). (3.23) 

(2) 

We thus see that the 2^1 contribution not only generates the "fluctuation" term, /c2— >i(^)) 
but it does also modify the (t^) term. If this term is to be unaffected by the additional 
vertex, as in (2.5), then we have to set k2_^i = 0. This is actually very similar to what 
happens in the full model. In that case the integral over the proposed 2 — > 1 vertex has 
to be zero, which implies that the vertex cannot be positive definite, as was noted in [19]. 
We can also try to add another vertex such that the total contribution to the (t^) term 
cancels. Assume for example the existence of an additional 2^0 vertex g2~>o- We then 
get 

drit') = (4'ii + + (/S2 - - &){t') - fi%{t') (3.24) 

from which we conclude that 

4'ii+5?io = (3-25) 
(2) (2) 

which means that either k2Si 5'2-»o negative. Thus we conclude that this 

(2) 

approach has big problems, as one must choose k2_^^ either to be 0, or it must be negative, 
which means that one cannot obtain a probabilistic formulation. 

3.4 The 1 + 1 dimensional toy model 

A somewhat more complicated 1+1 dimensional model is presented in [18]. The structure 
of this model is very similar to the 1+0 dimensional model, but the difference is that this 
time a dipole state is not only specified by the total number of dipoles, but it also depends 
on the distribution of these dipoles along some additional transverse axis. 

We denote the position of a dipole along this axis with Xi, and the generic n-dipole 
state is denoted = . . . , rE„). The assumption in [18] is that the dipole state 

evolves only by the addition of a single new dipole at some position Xn+i- In that case the 



frame independence equation in ( |3.5D can easily be solved, and the simplest solution for 



the total splitting rate TZi{{n}) is given by [18] 



-r,/r 1 r 1 , ^ 1 -n"=i(l -'^(x^lxn+i)) 

TZ{{xi} ^ {xi} + Xn+i) = — . (3.26) 
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Here r is a constant which can, by a redefinition of Y, set to be equal to 1. 

We now show that the evolution can once again be formulated as a probabilistic process 
in terms of coherent — > A; + 1 transitions as in sec |3.2| . In this case we assume we have a 
system of dipoles, Xi, which live on a one dimensional spatial axis. We assume this axis to 
represent the position of the "point-like" dipoles. This axis is assumed to be continuous, 
so that the index i actually represents a continuous label Xj. An isolated dipole Xi can 
then emit another dipole Xj at position xj with a probability nj = T{xi\xj). However, in 
the presence of more than one Xi, the new Xj can also be emitted coherently from several 
dipoles with a probability given by the product of the individual emission probabilities. 
For a system of n dipoles located at positions xi, . . . , x„, the total k ^ k + 1 splitting rates 
~^ i^i} + ^n+i) are then given by 

^ n n 
Tl'lli.^liixi} {Xi} + Xn+l) = -ry X] ^ii,n+in2,n+l ■ ■ ■ Ti^^,n+1 JJ ~ Tim,n+l), 



n^---^ifc m^l, 



(3.27) 



and their sum satisfies 

n n 
JZ^tlk+lii^^} ^ {^^} + ^n+l) = 1 - 11(1 - (3.28) 



k=l 



which is equal to ( |3.26| ). Note that once again the positive definite splitting rates do not 
only depend on the state of the emitting dipoles, but also on the state of the other dipoles 
in the cascade. This is unavoidable if one wants to obtain a probabilistic evolution. 

The A; — > /c + 1 splitting rates in ( 3.27| ) are very similar in structure to the processes 



generated by the dipole swing, to be discussed in the forthcoming sections. Anticipating 



the discussion there, we can interpret ( |3.27D as a process where the newly produced dipoles 
swing multiply with the rest of the dipoles in the cascade (a concrete example of this is 
shown in fig ^). Note that in the toy models both the splitting and the scatterings of the 
dipoles are determined by the quantities Tij . If we would assume that these also determine 
the swing probability, then the k ^ k + 1 splitting rates in ( 3.27] ) would describe processes 



where the newly produced dipole i swings with k — 1 dipoles from the cascade, and the 
factor rij(l ~ "^ij) could then be interpreted as the probability that i swings with no more 
than k — 1 dipoles. 

One difficulty is, however, that the swing in its form in the full model cannot really 
provide saturation in the toy models since the dipoles have no size here ^. This follows 
from the fact that both toy models have trivial topologies, in the sense that the dipole state 
is assumed to evolve only by the addition of a new dipole without changing the emitting 
state. In the toy models saturation occurs because k dipoles emit a single dipole with 
the same strength as a single dipole. In our implementation of the dipole swing however. 



''In [18], the spatial axis was interpreted as being related to the dipole size while we feel a more close 
analogy is to interpret it as a spatial coordinate where dipoles of some fixed size live. Irrespective of the 
interpretation, however, direct comparison with the full model is made difficult by the assumption that the 
toy model state only evolves by the addition of a single dipole. 
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saturation occurs because the swing decreases the sizes of the dipoles, and smaller dipoles 
have a smaller probability to interact. 



4. Evolution in the full model 

We will here argue that the evolution in the full model also can be formulated as a proba- 
bilistic process in terms of A; — > A; + 1 transitions. 
The S'-matrix is in the full model given by 

N M . X 

Sy = Y. Pn{Yo)Pm{Y -Yo)llll(l- /(x„ y^\xJ,yJ) J . 

N,M i=l j=l ^ ^ 

N M > 

(4-1) 

i=i j=i I 

where fij is given by (2.3). 

Let us consider the evolution initiated by a pair of oppositely moving qq pairs. We 
then consider generic dipole states containing — 1 gluons, which at each rapidity 
step can evolve into states containing N gluons. The splitting rate is denoted 

The scattering between the states £/n and ,'3^m is then frame independent if 

E /^(^M-^Am){i- n{i-( E - E 



(4.2) 



The notations in this equation are as follows. The integral J"^ denotes the integration over 
the transverse position of the A^th emitted gluon. In Yljes^m^ index j runs over all 
dipoles in the state ^m- The set denoted by -s/n+i \ s^N consists of those new dipoles 
produced in the last step of the evolution. Similarly, j^/at \ s^n+i denotes the set of all 
dipoles which are present in s^n, but not in i.e. those dipoles which were removed 

from the cascade in the last step of the evolution. Finally, f{i\j) stands for the scattering 
amplitude between the dipoles i and j, i.e. the expression fij in ( |2.3D and ( [4.11) . The sum 
X^^/jv+i ^ gluon states which can be reached from in one step. 

In the toy model analogy, the difference EiG^/iv+i\.c/jv ~ Ejei/jv\.£/jv+i consists of only 
the newly produced dipole, i, since all others dipoles are assumed to be unaffected by the 
evolution. In that case eq (^) reduces to 



(4.3) 
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The most simple solution is given by 

which we recognize from (1^. We will now see that the evolution in the full model can 
be formulated probabilistically in terms oi k ^ k -\- \ vertices as in the toy models. At 
first we will formulate the evolution as in eq ( |3.7D which implies that the 2^3 transition 
will appear to have negative sign. However, we know from above how to treat these signs, 
and thus give the evolution a probabilistic interpretation. We will also see how one can 
interpret these vertices in terms of the dipole swing. 

To this end, we consider first the situation where N = 2 and M = 1. The state 
must consist of two connected dipoles since we know that a single, isolated dipole (the state 

evolves by a dipole splitting (assuming the evolution is initiated by a single qq pair). 
We then denote the two dipoles in £^2 with a and b, and the single dipole in the state 
is denoted by /. We can then write ([4.2]) as 



^ fn{a,b^^3)i - Yl )fm = [ni^ih,z) + {z,i2))x 

X {^{l, a, z) + .^{l, b, z) - ^{l, a, z)^{l, b, z)^ (4.5) 

where 

^(/, a, z) = f{h,z\a) + f{z, l2\a) - f{l\a). (4.6) 

Here li and I2 denote the transverse positions of the partons of dipole I. Prom studying 
the case N = M = 1 (the scattering between two elementary dipoles), we know that each 
isolated dipole evolves by a 1 — > 2 splitting. We therefore write Yl^/^ 7^(a, b — =2/3) as 

Y'^ia,b^£/3) = n^^\a^ {ai,z) + {z,a2)) + n^^Hb^ ibi,z) + {z,b2)) 
^3 

+ 5:7^(^)(a,6^^f). (4.7) 



4'' 



(2) 

Here denotes the set of all 2 gluon states which can be reached from ^ = {a, 6} via 
the vertex TZ^'^\ Now, we know from [11] that the incoherent contributions, TZ^^\ above 
are equal to the first order contributions (in f{i\j)) in ([4.5|). Thus we are left with the 
equation 



n{l ^ {h,z) + (z, h)) (^^(/, a, z) ■ ^(/, 6, z)J . (4.8) 

In the toy model, where the dipole state evolves by the addition of a single dipole i only, 
we know that Ti^'^\a,b ■s^^'^^) = ~ f {^W) f {^W) (in case we use the formulation in (|3.71)). 
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Indeed in that case we see that both sides in (|4.8| ) equals — f {i\a) f {i\b) f . It is then 
clear that we must have a 2 ^ 3 transition. 

By similarly studying the case where = 3 and M = 1, we would conclude that we 
need an additional 3 — > 4 dipole vertex and so on. We are thus led to a picture where 
the dipole state evolves by /c — > A; + 1 transitions. If these transitions can be generated by 
combining the dipole splitting with the dipole swing, as we will argue in the next sections, 
we furthermore obtain a probabilistic interpretation of the evolution, as was discussed in 
the end of sec 

Before we going on, we also note that care has to be taken to the fact that the frame 
independence equation in ( |4.2| ) may contain divergences. In the original dipole model, 
these divergences arise from the dipole splitting kernel in but the frame independence 

equation is still finite since the expressions in the brackets in (^]^) vanish at these singular 
points. This is both due to the topology of the dipole splitting and also to the colour 
transparency of small dipoles. Any new vertex to be introduced into the model must 
retain this property since otherwise equation ( [4. 2]) would not make any sense. 



5. The Colour Topology of the Evolution 
5.1 Colour Flow 

Although we are not able to explicitely write down the splitting rate TZ^^^a, b ^^^^) in 
the full model, we will in this section argue that the correct topology of the evolution is 
the one induced by the dipole swing. 

We first write eq (|4.2|) in Mueller's original formulation. 



N . M . 

/ M{i\z)l 1 - 11(1 - ( E -Y^ )f{k\j)) = 

i=l ^ ^ j=l kenew k&old ^ 

M.N . 

= 5] i-n(i-( E -Y.)fm)\- (5.1) 

i=l ^ ^ j=l kenew kGold 

Here Ai{i\z) is the usual dipole kernel in ( ^TiP for a dipole i emitting a gluon at posi- 
tion z, and for simplicity we denote with ^^k^newioid) the sum over the dipoles produced 
(destroyed) in the last step. 

In case the newly produced dipoles only scatter against one target dipole, eq ( ^Tl] ) 
linearizes, and in that case the equality is known to hold [11]. Note that this case does not 
restrict us to one pomeron exchange, it is only the dipoles produced in the last step which 
should scatter against a single dipole. There may be still be several scatterings between 
dipoles produced earlier. In such events, the colour topology may be described by open 
chains stretching between the target and the projectile. An example is shown in fig |l|. Here, 
the simple 1 — > 2 splitting is sufficient for producing all possible colour configurations. 

In case the newly produced dipoles scatter off two or more target dipoles, however, 
more complicated topologies are formed, and the simple 1 — > 2 splitting is not sufficient 
anymore. Let us ffist consider the case where each dipole is restricted to single scattering 
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Figure 1: A process where the newly produced dipoles exchange only a single gluon with the 
target. In this case the colour correlations can be generated by the dipole splitting only. Here the 
position of z is integrated over and is therefore not fixed, and we only show one possible colour flow. 
Note that each interaction implies a change in the colour flow, which goes from colour to anticolour 
as indicated by the arrows. It is not important whether there are also other interactions or not. In 
fig (A) we assume an additional interaction between (2:2, 1/2) and {u2,V2)- In fig (B) we show the 
same colour fiow in the corresponding Feynman diagram. 



only. In this case, the newly produced two dipoles can at most scatter off two target dipoles. 
An example is shown in fig Q where three different colour configurations are formed (there 
are three more configurations which can be formed by reversing the colour flow in each 
line). Here two right-moving dipoles {xi,yi) and {x2,y2), are connected to a single left- 
moving dipole {u,v). With the restriction that each dipole scatters only once, the dipole 
(n, v) can obviously not be connected to more than two oppositely moving dipoles. In this 
case one can produce all colour configurations by combining the dipole splitting with only 
one swing. This is in agreement with the findings in [16]. In fig ^ we show how the three 
configurations in fig ^ can be generated from a dipole splitting and one swing, by first 
forming the configuration in fig |^(^). Figure ^ can be compared to eq ( |4.8| ) where the 
right hand side of that equation describes the evolution and the scattering of the dipole 
{u,v) {I in ([4.^)). The vertex TZ^'^^ on the left hand side would then correspond to the 
diagrams showing the evolution of (xi, yi) and (x2, 2/2) into the three dipoles (xi, z), {z, ^2) 
and (x2,yi). 

We conclude that in an approximation where multiple scatterings are allowed, but 
with the restriction that each dipole scatters only once, the maximal correlation induced 
between the dipoles is that between a pair, and such a correlation can be generated by a 
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simple swing. 

Actually, in this approximation, explicit frame independence in zero transverse di- 
mensions can be achieved by including a 2 ^ 1 vertex in addition to the usual 1^2 
splitting. This reflects the fact that the maximal correlation induced is that between a pair 
of dipolcs. Note also that, in the situation described above, it is always only one out of 
the three dipolcs produced via the combination of the dipole splitting and the swing which 
interacts with the target. In that sense the swing corresponds to a 2 ^ 1 transition. We 
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Figure 3: The three configurations marked by (a), (6) and (c) in fig || can all be generated when 
the dipole {u, v) interacts with one of three dipoles from the configuration shown in fig (A). This is 
illustrated in fig (B). As illustrated, the configuration in fig (A) can in turn be generated via the 
dipole swing. 

can thus obtain an effective 2^1 transition without actually decreasing the number of 
dipoles. 

If a single dipole can scatter multiply, one swing will not be enough. In this case one 
dipole can for example split into two new dipoles, and these two dipoles can then interact 
with more than two target dipoles, inducing higher order correlations. Before going on, we 
note that there is an ambiguity in the statement that one dipole scatters multiply. Since 
each scattering implies a recoupling of the colour flow, a dipole which interacts is replaced 
by a new dipole. What we rather mean here is that the partons of the dipole can exchange 
multiple gluons. 

Consider the diagrams shown in fig ^. Here a single left-moving dipole {u, v) is linked 
to three right-moving dipoles (xi, yi), (x2, 2/2) and (x3,y3), as shown in the figure. The 
two colour configurations shown in the figure can then be generated as illustrated in fig 
Hj: First one generates the configuration {(X3, 1/2)5 (x2, yi), (xi, z), (z, ys)} by combining a 
dipole splitting with two swings as shown in fig ^{A) (this is obviously not the only process 
from which this final configuration can be generated). One of the four dipoles in this 
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Figure 4: The colour flow in a process where three right-moving dipoles {xi,yi), {x2,y2) and 
(xa, j/3) interact with a single left-moving dipole {u,v) after the emission of a gluon located at z. 
For simplicity, we only show two, marked (a) and (6), out of the four possible configurations. 



state can then collide with the dipole {u,v). If for example {u,v) collides with (xi,z), the 
configuration marked by (a) in fig ^ is produced. 

The same process can also be viewed as an evolution of the dipole {u,v), which then 
splits into {u,z) and {z,v), and fig ^ shows also how the two configurations in fig ^ can 
be generated when the dipoles {u,z) and {z,v) interact with (xi, yi), (x2, y2) and (x3,y3) 
exchanging now 3 gluons (fig ^(C)). Thus at least one of the dipoles {u,z) and {z,v) 
must scatter multiply in this case, since it would otherwise be impossible to generate 
the necessary colour correlations. As remarked in [12], the evolution equations in sec 3.3 
actually describe such events where a newly produced dipole scatters multiply. We can 
also compare the processes in fig |5| to eq (|4.2| ) where one side of the equation describes the 
multiple scatterings of the dipoles (n, z) and {v,z), while the other side describes how the 
3-dipole system evolves into a 4-dipole system which then exchanges a single gluon with 
{u,v). 
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{C) 

Figure 5: The configurations marked by (a) and (6) in figure ^can be generated when the dipole 
(m, v) interacts with different dipoles from the configuration in fig (A), as illustrated in fig B. There 
are two more configurations which can be obtained when (u, v) interacts with the other two dipoles 
in fig (A). In fig (C) we show how the same configurations can be generated when the evolution is 
instead put into the dipole {u,v). 

5.2 Generating arbitrary correlations using at most — 1 swings 

We might then expect that the correlation induced by the scattering between a single 
left-moving dipole and k right-moving dipoles can be generated by a splitting followed 
by /c — 1 swings in the right-moving system. Note that such a process gives a k ^ k + 
1 transition. In the example above this was possible since we were able to form the 
configuration {{x3,y2), (x2, yi), (xi, z), (z, ys)} by combining one splitting with two swings. 
This will always be possible if, given an arbitrary set of N dipoles we always can generate 
all possible (A^-l-l)-dipole states, by combining a splitting with at most A^ — 1 dipole swings. 
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We will below argue that this is indeed the case. In case we have N spatially disconnected 
dipoles, the proof is easy. However, starting from a single qq pair, the evolution does not 
generate spatially disconnected dipoles, and in this case the result is a conjecture. 

5.2.1 Spatially disconnected dipoles 

As a warm up, we first show the statement in case we have N spatially disconnected dipoles 
{{xi, yi)}^i- We then wish to evolve this state into some arbitrary N + 1 dipole state, 

) yp{i )), (5.2) 

using at most N — 1 dipole swings. Here p{i) is a permutation oi i = 1, . . . , N. We first 
start by emitting gluon z from the dipole {xk,yk), 

Y{{xi,yi) {xk,z){z,yk)Yl{ Xi, yi ). (5.3) 

i i^k 

The result then follows if we can show that, for an arbitrary permutation p{i), we can with 
— 1 swings always make the transformation 

TV TV 

Ylixi, yi) Ylixi, yp(i)). (5.4) 
1=1 i=l 

To this end, we perform the following swings in the indicated order, 

1 ixN,yN)iXp(N),yp{N)) ixN,yp(N)){Xp(N),yN) 

2 {xN-l,yN-l){Xp(N-l),yp{N-l)) {xN-l,yp{N-l)){Xp(^N^i),yN-l) 

etc. Then, after at most N — 2 swings we are either finished, or we have 

TV 

{xi , yp{2) ) {x2 , ) n ' yp{i) ) • (^•^) 

i=3 

We then need only one more swing (xi, yp(2))(x2, 2/p(i)) (^ci, yp(i))(a;2, yp(2)), and so after 
at most N — 1 swings we are finished. 

The problem is that generally the dipoles are not spatially independent, and one then 
has to be careful in performing swings, since they might generate zero size dipoles, i.e. 
colour singlet gluons which cannot be allowed. 

5.2.2 Dipole states initiated by a qq dipole 

Representation of the dipole states and the swing in terms of permutations 

Consider the evolution initiated by a qq colour dipole. In the original formulation of 
the dipole model, the dipole state at each rapidity Y consists of an open chain, C, of colour 
dipoles which are linked together via the gluons. Note the dual role played by the gluons 
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(0 aja2...aj^_j) 







(«k.«k,+l«ki+2) 



Figure 6: A generic dipole state formed after a rapidity evolution of Y , starting from a qq pair. 
The initial quark and the antiquark are both denoted by while the gluons are denoted by as 



explained in the text and in eq (5.8). The arrows on the dipoles indicate the colour flow, which 



goes from colour to anti-colour as before. 



and the dipoles, each gluon links together two dipoles, and each dipole links together two 
gluons. 

The inclusion of the dipole swing generates closed dipole loops, C, in addition to the 
open chain, C. The swing induces the transformations 



C ^ Ci + £2- 



(5.6) 
(5.7) 



In what follows, we will denote each A^-dipole state as an element of the permutation group 
Vn- For simplicity we suppress the transverse coordinates in the notation, and each gluon 
is denoted by a number indicating the order in which it was emitted, the first emitted gluon 
is denoted 1, the second 2 and so on. The initial qq pair is simply denoted by 0. A generic 
A^-dipole state containing — 1 gluons, with fco — 1 gluons in the open chain, and the rest 
in m closed loops each containing ki gluons, is denoted 



£/n = (Oai . . . ako-i){ako ■ ■ ■ ako+ki-i) ■ ■ ■ (av-m-i , , . . . aN-i) G Vn- 



(5.8) 



Here {ai}^Si is a permutation ofi = l,...,A^ — 1. The generic dipole state in (5^) is 
illustrated in figure Each arrow indicates the colour flow, and in the group theoretical 
notation in (5.8), each gluon points to the gluon to the right of it. The open chain 
is always represented by the cycle containing the element (the qq pair), and each cycle 
in (5.8) corresponds to a colour singlet. Since we cannot have colour singlet gluons, the 



numbers Ui cannot appear as 1-cycles. The only 1-cycle allowed is (0), which corresponds 
to a dipole formed by the initial qq pair. 

Every element in the group Vn belongs to a certain class, which is determined by the 
cyclic structure of the element. The group Va has 5 classes: 1111, 211, 31, 22 and 4. Here 
each Ti-cycle is represented by the number n. The state -s^n in (|5.8D belongs to the class 
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Figure 7: Illustration of eq (5.9). 



ko ki . . . km- The identity element is the permutation which takes every number onto itself, 
and has the cyclic structure 11 ... 1. 

A swing operation can be represented by an element of Visf which consist of one 2-cycle 
and {N — 2) 1-cycles, i.e. by an element belonging to the class 211 ... 1. Thus for example, 
the swing illustrated in fig ^ is represented by S{ai,aj) = (ajOj) Ilfc^i j('^fc)) and we have 

S{ai, Uj) (g) (. . . Oi-iOi . . . Oj-iOj ...) = {aiOj) (ok) (. . . Oj-iaj . . . Oj-iaj • • • ) 

= (. . . Ui^iOj ...){ai... Oj-i). (5.9) 

Here ig) denotes the group multiplication. The action of S{ai,aj) makes Oi^i point at aj, 
and aj-i point at ai, leaving all other ak unchanged as shown in the figure. 

Due to the fact that not every swing leads to a physically acceptable state, the number 
of allowed swings for a state containing A'^ dipoles is not simply ^N{N — 1). This would 
e.g. be the number of pairs in a reaction-diffusion type of formalism. In most formulations 
this is not taken into account, but we here wish to emphasize the importance of keeping 
track of the correct topology of the evolution. While this is not important in the original 
formulation of the dipole model where the dipole state evolves through the 1^2 splitting 
only, it is very necessary for transitions involving more than one initial dipole. In the 
appendix we show that the number of physically possible states Md, and the number of 
possible swings Ms are for N dipoles given by 

Mo{N) = (AT _ 1)! ^ L_l(Ar _ z), (5.10) 
1=0 

Ms{N) = ^{N -1){N -2)+n2, (5.11) 
where n2 is the number of closed loops containing 2 dipoles. 

Multiple swings in the N ^ N + 1 evolution 

The classes of the group Vn are connected to each other via the swing as illustrated 
in fig where each line means that two elements from the respective classes can be trans- 
formed into one another using one swing. Note that the longest distance is that between 
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Figure 8: Class diagrams for P4 and P5. 




Figure 9: Class diagrams for the subset of physical states of and Ve- The open chain is marked 
by *. 

4 and 1111, which requires 3 swings. In V^, we need 4 swings to go from 11111 to 5, as is 
also shown in fig ^ 

Generally, for Vn, any element in the class can be reached from the identity element 
using — 1 swings. Explicitely, we can write the A^-cycle (ji . . -Jn) as 

(ii • ■■jn) = S{ji,jN) (S> S{ji,jN^i) (S>---(E) S{ji,j2) (5.12) 

This also implies that, given any arbitrary element a E Vn, we can reach any other element 
b G Vn using at most — 1 swings. This is so since we can always find — 1 swings such 
that their product equals ba~^. 

However, not all classes fall into the subset of physically acceptable states, which in 
particular does not contain the identity element. Therefore we cannot a priori say whether 
or not the result above also holds for this subset. In fig |9|, we show the class diagrams 
of physically acceptable states for and Vq. Here n* denotes the open chain containing 
n — 1 gluons. Thus using this notation we would say that ^/^ in ( ^. belongs to the "class" 
k^ki . . . km- With a slight abuse of nomenclature, we will for simplicity continue to refer 
to these quantities as "classes", even though they do not constitute classes in the group 
theoretical sense. 
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Figure 10: Class diagrams representing the evolutions = 4 ^ A'' = 5, and A'' = 5 — > iV = 6. 
Here the double lines represent the dipole splitting. Note that for these we can only go in one 
direction, from an un-circled class to a circled one. 

Actually, the dipole splitting can be represented by the same class of elements as the 
dipole swing. Assume we are in the state s^fq. We then regard the splitting as a two-step 
process; first, we add the A^th gluon as a 1-cycle into the state =2^, formally writing 
as an element of Pat+i, 

i24r = (Oai . . .afco-i)(«fco • • • «fco+fci-i) • • • ("y-'-ifc- • • • aAf-i)(A) e Vn^\- (5.13) 

We put a bar on since, written in this way, it is not a physically acceptable state. Then, 
in the second step, we represent the emission of A^ from the dipole spanned between ai and 
Oj+i by operating on with <S'(A, aj+i) (see eq ( |5.9D ), since in that case (. . . a, a^+i . . . ) 
is replaced by (. . . A Oj+i . . . ). 

The class diagrams for the generic A^ ^ A + 1 evolution can be drawn in a similar 
fashion as before. In fig[l^, we show examples for A^ = 4 and N = 5. Here only the circled 
classes are physically acceptable, and it is one of these that we must end up in, starting 
from one of the un-circled ones. The maximal distance between any two circled classes in 
an A ^ A^ -|- 1 evolution is 2 for A^ = 3, 3 for A^ = 4, and also 3 for N = 5. Thus in 
this case this distance is not equal to A — 1 for a A-dipole state. This does, however, not 
automatically imply that we can reach any given state in less than A — 1 swings. 

We here conjecture that one can also for the subset of physical states go from £/iy to 
any £/n+i by combining a dipole splitting with at most A^ — 1 swings. In the appendix we 
show explicitely that this statement is true for A = 4 (the cases A = 1, 2 are trivial, and 
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A'" = 3 can be checked very easily). We have also checked this result for = 5, = 6 and 
N = 7, but we will for simplicity not present the calculations for these cases. In a process 
where k — 1 swings take place, k dipoles in the cascade get replaced by A; + 1 dipoles, thus 
giving a — > + 1 transition. 

The only exception to the statement above is when there are states containing an 
isolated triangle. For example, if we wish to go from the state (0)(12 3)(4) in 1*31, to the 
state (04) (13 2) in 2*3, we need 5 steps totally, 1 splitting and 4 swings. Thus for these 
states, — 1 swings are not enough (one needs A^ + 1 swings). This is directly related to 
the fact that the physical set of states does not allow the steps: (12 3) (1)(2 3) ^ (1 3 2). 
Therefore we need to use 4 swings rather than only 2 swings, which implies that we generally 
need N + 1 swings for these states. Note, however, that this problem does not appear for 
higher order cycles. The 4-cyclc (12 34) can for example be transformed into (1432) 
easily: (12 3 4) ^ (14)(2 3) ^ (14 3 2). 

However, we also note that the only difference between the configurations (12 3) and 
(1 3 2) is in the orientation of the dipoles. Moreover, the states = ("i "2 C(3)^N-3 
and i^l^ = {ai as a2)^N-3 have exactly the same weights in the cascade evolution. They 
are therefore always produced equally, and it is therefore not a problem if we cannot go 
between them using A^ — 1 swings. Finally, we note that the semi-classical approximation 
represented by the cascade evolution cannot take into account all quantum-mechanical 
interference effects. The quantum-mechanical states corresponding to the configurations 
{aia2a3) and {aia3a2) have colour factors Tr(rT^r^) and Tr(r»r^r^) respectively, 
and for finite Nc these states are not orthogonal. Although the interference is suppressed 
by l/N^, it is enhanced in for example the decay process T 3g, and is in this case quite 
large [22], which is also confirmed experimentally. 

6. Conclusions 

Mueller's dipole model gives a simple picture of the small-a; evolution which is also very 
suitable to use in a MC simulation. While it is known that it gives the correct evolution 
for dilute systems, a fully consistent version for dense systems, where saturation effects 
during the evolution cannot be neglected, is not known. There have been some attempts 
to interpret these saturation effects in terms of dipole mergings but it has not been possible 
to present a consistent probabilistic formulation. 

A consequence of neglecting the saturation effects during the dipole evolution is that 
the model is not frame independent. In a previous paper [15] (see also [17] for a more 
detailed account) we demonstrated that approximate frame independence can be achieved 
by including a so called dipole swing in the evolution (the swing was also suggested in [16] 
as a mechanism to generate pomeron loops). Based on this, we constructed a phenomeno- 
logical model which, implemented in a MC simulation, gives an almost frame-independent 
formalism. 

It has been quite difficult to analytically derive the relevant dipole interactions which 
would give rise to saturation effects in the dipole model in a way consistent with boost 
invariance. A very simplified treatment of the dipole evolution is offered by the toy model 
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introduced in [3, 11], and later also studied in [12, 13]. In this model it is possible to 
modify the evolution so that the formalism is explicitely frame independent. The evolution 
proceeds here by the addition of a new dipole at each step, in such a way that the total 
splitting rate saturates as the dipole occupation number gets large. As discussed in [12], 
this is actually quite similar to the way saturation occurs in the CGC formalism. 

In this paper we have first shown that it is possible to give a probabilistic interpretation 
to the toy model evolutions in terms of positive definite — > A; + 1 transitions. These tran- 
sitions describe the coherent evolution of the dipoles, which is an unavoidable consequence 
of the requirement that the transition rates be positive definite. The evolution can also be 
formulated in a more close analogy with a standard reaction-diffusion picture, where the 
k ^ k + 1 transition rates only depend on the k dipoles involved in the transition. In this 
case, however, these rates appear with alternating signs which implies that a probabilistic 
treatment is not possible. 

In the real dipole model such positive definite vertices can be generated by combining 
the dipole splitting with the dipole swing. In a A: ^ A; -|- 1 transition, a splitting is combined 
with k — 1 simultaneous swings. In the approximation where each single dipole only scatters 
once, we have seen that it is enough to combine each splitting with a single swing in order 
to generate the necessary the colour correlations. 

When each single dipole is allowed to scatter multiply, one needs to include more 
than one simultaneous swing. In this case the evolution proceeds by the k ^ k + 1 
transitions as in the toy models mentioned above, and we have further shown that for a 
system of N dipoles, one needs at most — 1 simultaneous swings in order to generate 
all colour correlations induced by the multiple dipole interactions. We therefore obtain 
a close analogy with the toy model evolutions, and the dipole swing furthermore gives a 
probabilistic interpretation of the evolution. This is easy to show for spatially disconnected 
dipoles, but it is also the case in the more relevant situation when the dipoles are connected 
in chains. 

This statement is strictly speaking not true for states containing a triangular loop, 
where only one orientation of this loop can be reached using at most N — 1 swings. This is, 
however, no problem because the two possible orientations always appear with the same 
weight. 
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A. The Number of Dipole States 



In this section we will demonstrate that the number of possible states for a system con- 
taining N — 1 gluons, together with the initial qq pair, is given by formula (5.10). 



We start by considering n gluons in a closed topology, i.e. a state containing one or 
more closed dipole loops. If all n gluons are in the same loop we obviously have (n — 1)! 
possible states. Next we might have n gluons in two loops. The number of such states 
is given by the number of elements in Vn which consists of one k- and one (n — A;)-cycle. 
There are ^ (^) (fc — l)!(n — A: — 1)! such elements. The symmetry factor 1 /2 comes from the 
fact that we can write the A:-cycle either to the left or the right of the (n — /c)-cycle. For a 
closed topology consisting of m loops, each containing ki dipoles, the number of possible 
states is given by 

fJi V kj J f \ ml ni=i ki 

where 



'^ki = n. (A.2) 



For each fixed closed topology with n gluons we also have {N — 1 — n)! possible states in 
the open chain. To write down the total number of states it is convenient to introduce 
a generating function G{z) whose series expansion gives the dipole state multiplicity. We 
have 



//V1\ ^ 1 / °° k\m oo 

^ ^ m=0 ^k=2 ^ fco=0 



where k^ is the number of gluons in the open chain. We also demand that each closed 
loop contain at least 2 dipoles as we do not allow colour singlet gluons. The constraint 
X^i^o ki = N — 1 is automatically ensured since we are looking for the (A^ — l)th coefficient 
in the expansion of G. The expansion of G gives 

G{z) = (AT - 1)! ^ ^ i^(M - ^ + l)z^. (A.4) 

M=0 1=0 

We then immediately see that the M = N — 1 coefficient is equal to eq ( ^.lOf) . Notice 
also that for large N the number of states approaches This can be compared to the 
number of possible states for A'^ dipoles formed by A^ spatially independent charge-anti- 
charge pairs, which is A'^!, and to the number of states in a system consisting of a single 
open dipole chain, as in the original formulation of the dipole model, which is (A^ ~ 1)'- 



B. The Number of Possible Swings 

In this section we demonstrate that the number of possible swings for a system containing 
A^ — 1 gluons, together with the initial qq pair, is given by eq ( ^.11 ). 
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Assume again that we have m closed loops each containing ki dipoles (i = 1, . . . ,m) 
with ki ^ 2. The open chain contains /cq gluons, and thus Yll.=o ki = N — 1. Within each 
closed loop we then have 

^ m 

-^ki{ki-3)e{ki^3) (B.l) 

i=l 

swings. The theta function takes into account the fact that we cannot have any swings 
in a loop containing only two or three dipoles. The number of swings between the closed 
loops is given by 

-. m 

-J^kikj (B.2) 

since there are no restrictions in this case. The number of swings between the open chain 
and the closed loops is given by 

m 

Y,{ko + m, (B.3) 

i=l 

and finally, the number of swings within the open chain is given by 

ko{ko - 1) 



2 

The total number of swings is then given by 



(B.4) 



-\ Y,{ki{ki - 3) - ki{ki - 3)(5fc,2} + + 2(A;o + l)(iV - 1 - A;o) + k^ik^ - 1) \ 

= i|(iV - 1 - kof - 3(iV - 1 - ko) + 2n2 + 2{kQ + 1){N - 1 - fco) + A;o(A;o - 1)| 
= i(iV-l)(iV-2)+n2 (B.5) 
where n2 is the number of closed loops containing 2 dipoles. 

C. More Details on the ^ + 1 Evolution 

In this last appendix, we will explicitely prove that N — 1 swings are enough to reach any 
arbitrary state for N = A. We have also explicitely checked the cases N = 5, but we will 
for simplicity not present these calculations. We will very briefly try to sketch the case 
when N = 7. The cases N = 1,2 are trivial, and we also omit the case N = 3 which is 
very easy to work out. 
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C.l N = 4 

For this case, the class diagram is shown in fig |l^. Assume first that we are in the class 
4*1. An arbitrary element in this class is given by (Opi p2P3)(4), for some permutation 
{p{i)}- We must reach any arbitrary element using at most 3 swings, and thus using at 
most 4 steps, counting the splitting as one step. 

4*1 —» 1*22: We see from fig that we have to reach any element in 1*22 using at 
most 3 steps (or else we would need at least 5 steps). An arbitrary element in 1*22 can 
be written (0)(p'i P2)(P3?'4) where {p'{i)} is some other permutation. Without any loss of 
generality we might as well assume p'^ = 4. Then we can always start by putting 4 next 
to P3 in the step 4*1 ^ 5*. In the next step we can then always isolate {p'^4:) = {p'^p'^) 
in a 2-cycle. We then have an element (07r(l)7r(2))(p3 P4) where (7r(l), 7r(2)) = (^'1,^2) 
{P2,p'i). Finally we can separate(07r(l)7r(2)) (0)(7r(l)7r(2)) to reach {0){p'iP2){p'^p'4). To 
summarize, we can go through the following steps 

4*1^5*^3*2^1*22, (C.l) 

and reach any element in 1*22 using at most 3 steps. 

4*1 —> 3*2: Here we can go in either 2 or 4 steps. We then want to go to an element 
(Op']^ P2)(P3 ^4)- Again we start by putting 4 to the right of the number which appears to 
the left of it in (Op'^ P2)(P3 P4)- If this number is (i.e. iip[ = 4), we can go from 5* to 2*3, 
putting and 4 in 2*. Then we just go back to 5* to obtain {Op'i P2 p'^ p'^) or {0piP2p'4^p'^). 
In either we case we can split this chain into (Op'^ P2)(P3 ^4)- Thus we go through 

4*1^5*^2*3^5*^3*2. (C.2) 

If p'^ 7^ 4, we go from 5* to 3*2, isolating 4 and its partner in 2. If the three elements in 3* 
are not in the right order, we can split 3* into 1*2 and then go back to 3*. Thus we can 
go through the steps 

4*1^5*^3*2^1*23^3*2. (C.3) 

4*1 — *• 2*3: We here want to go to the element i0pi)iP2P3P'4)- Again we start by 
putting 4 together with its final partner (4 is always put to the right of its partner) . If its 
partner is 0, we can first go from 5* to 3*2 and then go back to 5* so that the elements 
(P2P3P4) have the correct permutation. Then we can in one step go the final configuration, 
and thus complete the process 

4*1 ^ 5* ^ 3*2 ^ 5* ^ 2*3. (C.4) 

If 4 is in the 3-cycle in the final element, we can go from 5* to 3*2 and then to 5* again 
to put all the elements together in the correct positions. This is possible since all three 
elements in the final 3-cycle must have the correct permutation. Then we can finish by 
going from 5* to 2*3. Thus we can choose the path 

4*1 ^ 5* ^ 3*2 ^ 5* ^ 2*3. (C.5) 

Here the 2-cycle in 3*2 contains 4 and its partner. 
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4*1 — > 1*4: Here we want to reach an element (0)(y^ P2P3P4)- Since we always put 4 
(without loss of generality we can assume p'^ = 4) next to its partner in the first step, all 
we need to do is to isolate them (pg and ^4) in a 2-cycle by going from 5* to 3*2. Then we 
go to 1*22, after which we can simply join 22 to 4, to obtain any desired state. We thus 
have the path 

4*1 ^ 5* ^ 3*2 ^ 1*22 ^ 1*4. (C.6) 

4*1 —s- 5*: We can here use at most 3 steps. The final element we want to reach has 
the form {^'p'xv'^v'z'P'a)- There are two cases, either 4 and its partner are linked to 0, or 
they are not. If they are, we can split 5* into 3*2 where 3* contains 0, 4 and its partner. 
Then the other two elements can always be put back in 5* in the right position, so that we 
reach any 5* element by 

4*1 ^ 5* ^ 3*2 ^ 5*. (C.7) 

In the second case, 4 and its partner are not linked to (they are P2 and pg). Then we can 
split 5* into 2*3 where 3 contains 4, its partner and one of the other two elements. They 
will automatically have the correct permutation, and we can then get the desired state by 
joining 3 and 2* into 5*. Then we have used the path 

4*1 ^ 5* ^ 2*3 ^ 5*. (C.8) 

We have thus seen that we can reach any arbitrary state in AT = 5 from 4* 1 by combining 
at most 3 swings with a splitting. 

Below wc list the cases where we start from 2*21. In this case we have an initial 
element (0pi)(p2P3)(4). By using a splitting first, we can either go to 3*2, or to 2*3. 

2*21 — » 5*: First we fix 4 and its partner as usual. If the partner is pi, we can separate 
(Opi 4) into (0)(pi 4), and then we can join (pi 4) with (P2P3) to obtain an element in 1*4. 
Then in one step we can go to the desired 5* state. If its partner is 0, and the other three 
elements do not have the correct permutation, we can isolate two of them in a 2-cycle (after 
putting 3* and 2 into 5*), and then put them back into the 5* state in the correct position. 
Thus we can through the two paths 

2*21^3*2^1*22^1*4^5*, (C.9) 
2*21-^3*2^5*^3*2^5*. (C.IO) 

If the final partner of 4 is either p2 or ^3, we first go to 2*3. Then it is easily seen that the 
two paths, 

2*21 ^ 2*3 ^ 5* ^ 3*2 5* (C.ll) 
2*21 ^ 2*3 ^ 5* ^ 2*3 ^ 5*, (C.12) 

can take us to any arbitrary element in 5*. 

2*21 — *• 1*4: If 4 is next to either p2 or p3, we can directly from 2*3 go to 5*, and then 
to 1*4. If 4 is next to pi, we can from 3*2 go to 1*22 and then in one more step we can 
reach any 1*4 state. Thus we can follow the paths 

2*21-^2*3^5*^1*4 (C.13) 
2*21 -> 3*2 1*22 1*4, (C.14) 



-29- 



to reach any state in 1*4 in maximum 3 steps. 

2*21 1*22: If 4 is next to pi we can finish in 2 steps, 2*21 3*2 -> 1*22. If 4 is 
next to either p2 or ps, we can first join 2*3 into 5*, and then isolate 4 and its partner in 
the 2-cycle in 3*2. Then we need only one more step. We thus have the steps 

2*21 ^ 2*3 ^ 5* ^ 3*2 ^ 1*22. (C.15) 

2*21 3*2: If pi is partner to 4, we need at most go to 1*22 from 3*2, and then back 
to 3*2 to finish. If 4 is next to 0, we can first go to 5* from 3*2, and then spht 5* into the 
desired 3*2 state. If 4 is next to p2 or ps, we can again go to 5* and then directly to 3*2. 
Thus we have the steps 

2*21 ^ 3*2 ^ 1*22 ^ 3*2 (C.16) 
2*21 ^ 3*2 ^ 5* ^ 3*2 (C.17) 
2*21 ^ 2*3 ^ 5* ^ 3*2. (C.18) 

2*21 — ^ 2*3: If 4 is in the final 2-cycle {i.e. next to 0), all we need to do is to join 3*2 
into 5*, after which we can extract the final 3-cycle in one step. If one the other hand 4 
is in the final 3-cycle, we are after one step either finished, or we can from 2*3 go to 5*, 
putting and its final partner together, after which wc can split 5* split into 2*3, obtaining 
the desired state. Thus we can go choose one of the paths, 

2*21 ^ 3*2 ^ 5* ^ 2*3 (C.19) 
2*21 ^ 2*3 ^ 5* ^ 2*3. (C.20) 

Finally, wc check the case when we start from 1*31. 
1*31 1*22: Here we only need two steps: 

1*31^1*4^1*22 (C.21) 

which can be easily seen. 

1*31 — » 1*4: This case is almost trivial, and we can see that we need at most three 
steps: 

1*31 1*4 ^ 1*22 ^ 1*4. (C.22) 

1*31 3*2: If 4 appears in the 2-cycle in 3*2, all we need is to take the steps 

1*31 ^ 1*4 ^ 1*22 ^ 3*2. (C.23) 

In the second step we here isolate 4 and its partner in one of the 2-cycles. If 4 instead 
appears in 3*, we can go through either 1*31 1*4 ^ 5* ^ 3*2, or 1*31 ^ 2*3 ^ 5* ^ 
3*2, depending whether or not 4 appears next to in the final configuration. 

1*31 5*: Here we can again have 4 either to the right of or not. If not, we just 
take the steps 1*31 1*4 1*22 — > 3*2 — > 5*. If it is next to 0, we instead take the steps 
1*31 2*3 ^ 5* ^ 3*2 5*. 

1*31 2*3: If 4 appears in 3, we can just go through 1*31 -> 1*4 ^ 5* ^ 2*3 and 
finish. However, if 4 is next to 0, A'^ — 1 swings are not enough to go to (04)(j)ip3P2) as 
we have already discussed in the main text. We have also already noted that this is not a 
problem for the frame independence. With this remark we finish the case = 4. 
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Figure 11: Class diagrams for the evolution N = 7 
lines between the un-circled classes. 



N = 8. Here for simplicity we do not draw 



C.2 iV = 7 



The class diagram for the evolution : 7 ^ 8 is shown in fig |11|. Here it is obviously too 
tedious to explicitely check all possible connections. However, one can now use the results 
from the previous cases to simplify the analysis. For example, let us consider the case 
where we want to go from 7*1 to 1*322. If 7 appears in the final 2-cycle, we can pull it out 
together with its partner after 2 steps, and the question is then whether we can go from 
6* to 1*32 in 5 steps, and we know from the case N = 5 that this is indeed possible. If 7 
appears in the final 3-cycIe, we can after 4 steps isolate it with its partners, and then we 
need to go from 5* to 1*22 in 3 steps, which we also know is possible. Another example is 
if we want to go to a final configuration in 2*42, where 7 appears in 2. After 2 steps, 7 and 
its partner can again be isolated, and we then need to go from 6* to 2*4 in 5 steps. Again 
we know that this is indeed possible. We can similarly work out the rest of the cases. It 
is also interesting to note that we can here go from 1*331 to 2*33 in 7 steps, even if the 
initial and final states differ by the orientations of the two triangular loops. 
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